Ultracold Superstrings in atomic Boson-Fermion mixtures 
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We propose a setup with ultracold atomic gases that can be used to make a nonrelativistic super- 
string in four spacetime dimensions. In particular, we consider for the creation of the superstring a 
fermionic atomic gas that is trapped in the core of a vortex in a Bose-Einstein condensate. We explain 
the required tuning of experimental parameters to achieve supersymmetry between the fermionic 
atoms and the bosonic modes describing the oscillations in the vortex position. Furthermore, we 
discuss the experimental consequences of supersymmetry. 
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In recent years, three topics have in particular at- 
tracted a lot of attention in the area of ultracold atomic 
gases. These topics are vortic es IH 0, HI d Boson- 
Fermion mixtures [fl 0, IS El llfij- and optical lattices 
[llL fl2j |. In this Letter we propose to combine these 
three topics to engineer a superstring in the laboratory, 
i.e., a line- like quantum object with both bosonic and 
fermionic excitations and a supcrsymmetric hamiltonian 
that is invariant under interchanges of these excitations. 
The physics of a vortex line in a one-dimensional opti- 
cal lattice has been studied recently [l3l llij . Because of 
the optical lattice, the transverse quantum fluctuations of 
the vortex line are greatly enhanced in this configuration. 
The vortex can therefore be viewed as a quantum me- 
chanical string and it forms the bosonic part of our super- 
string. In addition, we propose to trap fermionic atoms 
in the vortex core, to make a nonrelativistic version of 
a so-called Green-Schwarz superstring in four spacetime 
dimensions. Because our ultracold superstring is non- 
relativistic, it is not constrained to the ten-dimensional 
spacetimes in which superstrings are usually studied in 
high-energy physics. The precise mathematical connec- 
tion with string theory is currently under investigation. 

Apart from the connection to string theory, our ultra- 
cold superstring is also of interest in its own right. To the 
best of our knowledge, it is the first condensed-matter 
system proposed, where supersymmetry can be stud- 
ied experimentally. The intruiging possibility of observ- 
ing effects of supersymmetry is common in high-energy 
physics, but novel in a condensed-matter setting. In this 
particular case, the supersymmetry protects the super- 
string from spiraling out of the gas. This can be under- 
stood from the fact that the dissipation resulting in this 
motion has two sources, namely the creation of two addi- 
tional bosons in the transverse oscillations of the vortex, 
and the production of an additional particle-hole pair of 
fermions. At the supersymmetric point these two con- 
tributions interfere destructively and the stability of the 
superstring is greatly enhanced. Moreover, the ultracold 
superstring allows for the study of a quantum phase tran- 
sition that spontaneously breaks supersymmetry. Exper- 



imentally this will be directly visible by observing the 
superstring spiraling out of the center of the gas. Note 
that supersymmetry can only be realized at zero temper- 
ature and we consider only this case from now on. 

To form the superstring we start with a cigar-shaped 
Bose-Einstein condensate. The symmetry axis will be 
called the z-axis from now on. Rotation of the conden- 
sate along the z-axis creates a vortex. Above a critical 
external rotation frequency f2 c a vortex in the center of 
the condensate is stable. For fl < fl c the vortex is un- 
stable, but because of its Euler dynamics, it takes a rela- 
tively long time before it spirals out of the gas f^. flfll flfi| . 
We analyze in detail the case of f2 = 0, i.e., the situation 
in which the condensate is no longer rotated externally 
after a vortex is created. However, the physics is very 
similar for all Q < Cl c , where supersymmetry is possible. 

Next a one-dimensional optical lattice is imposed. 
Such an optical lattice consists of two identical counter- 
propagating laser beams and provides a periodic poten- 
tial for atoms. When applied along the z-axis of the con- 
densate, the optical lattice divides the condensate into 
weakly-coupled pancake-shaped condensates, each con- 
taining typically of the order of Nb — 10 3 bosonic atoms. 
Moreover, in the case of a red-detuned lattice, the gaus- 
sian profile of the laser beam provides also the desired 
trapping in the radial direction. In the one-dimensional 
optical lattice the vortex line becomes a chain of so-called 
pancake vortices. Since the quantum fluctuations of the 
vortex position are proportional to 1/Nb jl."^ . they are 
greatly enhanced in this configuration as compared to 
the bulk situation. An added advantage of the stacked- 
pancake configuration is the particle-like dispersion of 
the vortex oscillations, which ultimately allows for su- 
persymmetry with the fermionic atoms in the mixture. 
Nearest-neighbor pancake vortices attract each other due 
to the Josephson effect that is a result of the hopping of 
atoms between neighboring wells of the optical lattice. 
The stiffness of the vortex line is therefore determined 
by the hopping amplitude. The vortex oscillations are 
bosonic excitations with a tight-binding dispersion 
These excitations are known as Kelvin modes and pro- 
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FIG. 1: (color online) (a) Artist's impression of the setup. 
The disks represent the bosonic condensate density and the 
blue balls in the vortex core represent the fermionic density. 
The black line is a guide to the eye to see the wiggling of the 
vortex line that corresponds to a Kelvin mode, (b) Schematic 
of the setup. Here r is the radial distance in the xy-p\&ne. 
The pink and blue blobs represent the bosonic and fermionic 
densities, respectively. Moreover, A is the wavelength of the 
laser. The blue and red lines indicate the strength of the 
optical potential, respectively, for the bosons and fermions 
as a function of the z-coordinate. (c) Schematic fine struc- 
ture level scheme of the bosonic and fermionic atomic species. 
Because we consider only sufficiently large detunings the hy- 
perfine level structure is not resolved. 



vide the bosonic part of the superstring. Without the 
optical lattice Kelvin modes have already been observed 
experimentally [billisj]. 

For the fermionic degrees of freedom of the superstring, 
we need to place a fermionic atomic species in the vor- 
tex core. If the interspecies interaction between fermionic 
and bosonic atoms is repulsive, the fermions can be made 
to populate bound states in the vortex core where the bo- 
son density is lowest. More precisely, the ratio cibf/clbb 
of the boson-fermion scattering length and the boson- 
boson scattering length needs to be sufficiently large. 
The trapping of atoms in bound states within a vortex 
core has already been achieved for atoms in a different 
hyperfine state Q, 0] . The presence of fermions within 
the vortex core in principle provides a mass to the vor- 
tex. This modifies the Euler dynamics, but this effect 
can safely be ignored under all conditions of interest to 
us, because the ratio Nf/Nb of the number of fermionic 
atoms and the number of bosonic atoms per site turns out 
to be very small. The general properties of our proposed 
system are summarized in Fig. 1. 

A convenient choice for the boson-fermion pair is 87 Rb 
and 40 K, since such Boson-Fermion mixtures have re- 
cently been realized in the laboratory @, E El and 
because the resonance lines in these two atomic species 
lie very nearby. In this case, for typical bosonic den- 
sities, the ratio clb f / &b b needs to be larger than 2 to 



have a bound state in the vortex core. This is calculated 
by using the condensate density profile as an effective 
potential for the fermions. The interspecies scattering 
length can be made positive and large either by choosing 
the appropriate spin states or by means of various broad 
Feshbach resonances that can make the interaction re- 
pulsive while keeping the probability to create molecules 
negligible [l^. The calculations presented in this Letter 
are for the 87 Rb - 40 K mixture, but in principle it is also 
possible to use other mixtures. Another Boson-Fcrmion 
mixture that has been realized in the laboratory consists 
of 23 Na and 6 Li atoms .10] . This mixture is less conve- 
nient because the resonance lines are widely separated, 
so that the two species feel very different optical poten- 
tials and it is hard to trap both with a single laser. In 
addition, 6 Li is relatively hard to trap in an optical lat- 
tice because of its small mass. For these reasons, the 
23 Na- 6 Li mixture can only be used in a very restricted 
parameter regime, as shown later on in Fig. 2. 

The laser intensity needs to be strong enough that 
there is a bound state for the atoms at each site, but 
not strong enough to drive the system into a Mott insu- 
lator state [ill Il2l [23 , in which case the pancake vor- 
tices no longer bind together. Given the frequencies of 
the D\ and D2 resonance lines ojd x and lod 2 , respec- 
tively, the optical potential for the atoms is given by 
Vb,f{z) = Vb,f cos 2 [27rz/ A] , where the well depths obey 
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cu = 2nc/X is the laser frequency, and Qb.f are the Rabi 
frequencies for the bosonic and fermionic atoms in the 
mixture, respectively. 

Because of the use of a single laser the fermionic 
and bosonic atoms experience lattice potentials with 
the same periodicity, but with different well depths (cf. 
Fig. lb). Ignoring their interactions, the fermionic 
atoms and the kelvon excitations now both have a tight- 
binding dispersion 2J[1 — cos(feA/2)], with hopping am- 
plitudes Jf and Jk, respectively. Therefore, by mak- 
ing the hopping amplitudes of the kelvons and the 
fermionic atoms equal, the lattice laser can be tuned 
in such a way that the kelvon and fermion dispersions 
coincide, which is a first requirement for supersymme- 
try. For a sufficiently deep potential the atomic hop- 
ping amplitudes are given by the exact result Jb,f = 

HV b ,f) 3/4 {E b ,f) 1/4 cxp [-2^/V b ,f/E b ,f] /y/Z, where 

Eb.f = 2ir 2 H 2 /itib^X 2 are the recoil energies and tub^f 
are the atomic masses. The kelvon hopping amplitude 
is given by Jk = r[0, (£/R) 4 } Jb, where £ is the radial 
harmonic-oscillator length for the bosonic atoms, R is 
the Thomas-Fermi radius of every pancake, and r[0, z] 
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FIG. 2: Tuning of the lattice laser to obtain supersymmetry 
for 1000 bosonic atoms per lattice site: Rabi frequency for 
the bosonic atoms versus wavelength for 87 Rb- 40 K (full line) 
and 23 Na- 6 Li (dotted line). Note that for the blue-detuned 
part, i.e., A < 760 nm for the 87 Rb- 40 K mixture extra radial 
trapping is needed, either magnetically, or by using an extra 
running laser as discussed later in the text and shown in Fig. 
3. In Fig. 3 we display how to tune the running laser to 
obtain also supersymmetric interactions. The numbers (l)-(4) 
indicate the parameters for which this is calculated. The inset 
shows the average number of fermions per lattice site. This 
depends linearly on the ratio of the harmonic lengths in the 
axial and radial directions. This ratio should be sufficiently 
small to be radially in the Thomas-Fermi limit. For this plot 
a ratio of 1/5 is chosen. 

is the incomplete gamma function [l^. The ratio £/R 
is determined by the number of bosonic atoms per site 
Nb and the strength of their repulsive interaction. One 
criterion for supersymmetry is, therefore, 

T[0,(e/R) 4 ]J B = J F = t. (1) 

In Fig. 2 is indicated how to tune the wavelength and 
Rabi frequency of the lattice laser to achieve this. Note 
that we implicitly assumed that the fermionic hopping 
parameter is independent of the presence of kelvons, 
which is justified for the long wavelengths of interest, 
because then the difference between the positions of the 
nearest-neighbor pancake vortices is negligible with re- 
spect to the core size. 

A second requirement is that the chemical potential 
of the kelvons is equal to the chemical potential of the 
fermions. Therefore, the filling factor of the fermions 
has to be adjusted accordingly. Equating the chemical 
potential of the kelvons and fermions we get 

MA- = /iF = ^^(l-r[0,(£/i?) 4 ])^ M . (2) 

From this we can derive that the average num- 
ber of fermions per lattice site is given by Np = 
2 arcsin[-\/ /ip /4Jf]/tt. For typical parameters the num- 
ber of fermions per lattice site turns out to be of the order 



FIG. 3: Tuning of the additional laser to obtain supersymmet- 
ric interactions for 1000 bosonic atoms per lattice site. The 
Rabi frequency versus the wavelength of the running laser for 
different wavelengths of the lattice laser beam: (1) 400 nm, 
(2) 740 nm, (3) 850 nm, and (4) 1200 nm. See Fig. 2 for the 
corresponding Rabi frequencies of the lattice laser beam. 

of 0.1, as shown in the inset of Fig. 2. Tuning the sys- 
tem parameters to lie on the curves in Fig. 2 makes the 
string supersymmetric. The system is then invariant un- 
der unitary transformations of the bosonic and fermionic 
excitations among each other. 

In our supcrstring realization there are also boson- 
boson and boson-fermion interactions. The kelvons inter- 
act repulsively among each other when fl < fl c 0|. In 
addition, a repulsive interaction between the kelvons and 
the fermionic atoms is generated by the fact that physi- 
cally the presence of a kelvon means that the vortex core 
moves off center, together with the fermions trapped in it. 
Because of the confinement experienced by the trapped 
fermions in the radial direction, this increases the energy 
of the vortex. The kelvon-kelvon interaction coefficient 
isgiven by V KK = (T[0, {i/Rf] - ^)h 2 /2N B m B R 2 = U 
[14| . and the kelvon- fermion interaction coefficient is 
found to be Vkf = CfR 2 /2Nb- Setting these coeffi- 
cients equal to each other gives the third condition for 
supersymmetry, i.e., 

g = (i) 4 (n».w*)Vf). ,3, 

where Cb,f are the spring constants associated with the 
radial confinement of the bosonic and fermionic atoms, 
respectively. 

For the 87 Rb- 40 K mixture, supersymmetric interac- 
tions require the ratio of the radial trapping potentials for 
the bosonic and fermionic atoms to be large, of the order 
of 100. This cannot be achieved by a magnetic potential, 
since the magnetic moment of the two species is almost 
the same. To overcome this problem, we propose to apply 
an extra running laser beam along the z-dircction. The 
new laser beam docs not influence the one-dimensional 
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FIG. 4: The second-order Feynman diagrams describing the 
dissipation processes due to the creation of (a) additional 
pairs of bosonic and (b) fermionic excitations. The full lines 
correspond to bosons and the dotted lines to fermions. Be- 
cause the right diagram contains a fermionic loop, it comes 
with an additional minus sign and the two diagrams thus can- 
cel exactly in the supersymmetric case. 

potential wells, but it docs change the radial confinement 
and can therefore be used to tune the interaction terms 
equal. In principle this second laser also introduces inter- 
ference terms, but they turn out to be negligible, since 
they oscillate faster than the atoms are able to follow. 
Hence, the intensities of the two lasers can simply be 
added. The detuning of the second laser has to be op- 
posite to the detuning of the first laser. For the second 
laser, we can again independently choose both the wave- 
length and the Rabi frequency as shown in Fig. 3. This 
can be used to minimize the atom loss due to the red- 
detuned laser, but it turns out that atom loss is always 
quite small anyway for reasonable system parameters: at 
worst the lifetime of the system is already 3 seconds. 

Combining all these effects, our superstring is de- 
scribed by the supersymmetric hamiltonian 

H = -tYM*i + 4<h) (4) 

m 

+Y^(-iM b i - » c h + 2 b i b i b * b * + Ub i b i4 c i)- 

i 

Here bi is the annihilation operator of a kelvon at site i, 
Ci is the annihilation operator of a fermion at site i, and 
(ij) means that the summation runs over neighboring 
sites. Note that supersymmetry means that [Q,W] = 0, 
with Q = b i c \- For fi < O c , the chemical potential is 
positive and the kelvons are thus unstable towards Bosc- 
Einstein condensation, which physically corresponds to 
the tendency of the superstring to spiral out of the cen- 
ter of the gas. The metastable superstring at the center 
can be regarded as a system that is quenched deep into 
an ordered phase, but is yet to move to equilibrium. The 
Bose-Einstein condensation of kelvons corresponds to a 
quantum phase transition that spontaneously breaks su- 
persymmetry, which manifests itself by the kelvon dis- 
persion becoming linear at long wavelengths. Thus our 
proposed system naturally allows for a detailed study of 
the noncquilibrium dynamics of supersymmetry breaking 
by monitoring the dynamics of the superstring. 

It is important to realize that if supersymmetry is un- 
broken, i.e., the superstring is in a supersymmetric quan- 
tum state, the kelvon and fermion modes should have the 
same average occupation number, i.e., (£>J&») = (cjcj). 
Since the total number of fermions is conserved, the su- 
persymmetry in this case protects the superstring against 



dissipation, because the spiraling out of the superstring 
would imply the creation of extra kelvons. This results 
in a very long lifetime of the ultracold superstring, which 
can also be understood theoretically as an exact can- 
cellation of the second-order Feynman diagrams describ- 
ing the dissipation processes due to the creation of ad- 
ditional pairs of bosonic and fermionic excitations, re- 
spectively, which are shown in Fig. 0] In addition, 
the equality of the average kelvon and fermion num- 
ber allows us to devise an experimental measure for the 
proximity to the supersymmetric point. The number 
of kelvons {b\bi) can be obtained from the number of 
bosonic atoms Nb and the mean-square displacement 
(r ) of the pancake vortices, which can be measured by 
imaging along the z-direction the size of the circle within 
which the vortex positions are concentrated Ql- The 
average fermion number (eje,-) = Np can be determined 
by absorption measurements. Hence also the quantity 
((blbi) - (c\ Cl )) 2 = (N B (r 2 )/R 2 - 1/2 - N F ) 2 can be 
measured. This quantity has an absolute minimum of 
zero at the supersymmetric point, so that its magnitude 
is a measure of the deviation from supersymmetry. The 
experimental precision that can be reached in approach- 
ing the supersymmetric point will mainly be limited by 
the uncertainty in the total number of fermions in the 
system. Nevertheless, the observation of the effects of 
supersymmetry in our setup should be in experimental 
reach with the existing technology. 

We are grateful for helpful discussions with Randy 
Hulet and Jan Ambj0rn. 
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